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O-i. Abstract 



This paper presents the first analysis of diffractive photon dissociation events in deep 
inelastic positron-proton scattering at HERA in which the proton in the final state is 
detected and its momentum measured. The events are selected by requiring a scattered 
' proton in the ZEUS leading proton spectrometer (LPS) with xl > 0.97, where xl is the 

fraction of the incoming proton beam momentum carried by the scattered proton. The use 
of the LPS significantly reduces the contamination from events with diffractive dissociation 
of the proton into low mass states and allows a direct measurement of t, the square of 
the four-momentum exchanged at the proton vertex. The dependence of the cross section 
on t is measured in the interval 0.073 < |t| < 0.4 GeV^ and is found to be described 
by an exponential shape with the slope parameter b = 7.2 ± l.l(stat.)^Q'g(syst.) GeV~^. 

The diffractive structure function F^^"^^ is presented as a function of ~ 1 — xl and /?, 
the momentum fraction of the struck quark with respect to Xp, and averaged over the t 
interval 0.073 < |t| < 0.4 GeV^ and the photon virtuality range 5 < < 20 GeV^. In the 
kinematic range 4x 10~^ < x^ < 0.03 and 0.015 < f3 < 0.5, the Xp dependence oi F^^'^^ is 
fitted with a form (l/x^ , yielding a = 1.00 it 0.09 (stat.)lQ'Q5(syst.). Upon integration 

over t, the structure function F^^^^ is determined in a kinematic range extending to higher 
Xjp and lower (3 compared to our previous analysis; the results are discussed within the 
framework of Regge theory. 
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1 Introduction 



A distinct class of events in deep inelastic neutral current positron-proton scattering (DIS) is 
experimentally characterized by the proton emerging intact (or excited into a low-mass state) 
and well separated in rapidity from the state produced by the dissociation of the virtual photon 
[|l], 1^. These large rapidity gap events can be interpreted as being due to diffractive interactions 
mediated by the exchange of a colourless object with the vacuum quantum numbers, generically 
called the pomeron. 

In diffractive single dissociation DIS events at HERA, ep eXp, the virtual photon dissociates 
into a hadronic system of mass Mx, while the proton remains intact; for small |t| GeV^), 
where t is the square of the four-momentum transfer at the proton vertex, the scattered proton 
escapes through the beam pipe without being detected in the central apparatus. Therefore 
diffractive interactions have been studied so far in the HI and ZEUS experiments either by 
requiring a large rapidity gap 0, H or exploiting the different behaviour of the ln(M|^) 
distribution for diffractive and non-diffractive events In both approaches, in the ZEUS 
results, the hadronic mass Mx was limited to values below 20 GeV. 

This paper presents the first analysis of diffractive DIS events in which the proton in the 
final state is detected and its momentum is measured. The measurement was performed with 
the ZEUS detector at HERA, using the leading proton spectrometer (LPS) which detects the 
scattered proton at small angles (^1 mrad) with respect to the incoming proton beam and 
measures its momentum. The events were selected by requiring that the scattered proton 
carry a fraction of the incident proton beam momentum, xl, greater than 0.97, a region where 
pomeron exchange dominates. The use of the LPS allows a direct measurement of t, and 
extends the measurement to higher values of Mx than in our previous analyses. 



2 The diffractive structure function F^^^^ 

The dependence of the total hadron-hadron and photon-hadron cross sections on the centre 
of mass (cm.) energy is related in the Regge approach to the pomeron trajectory. A fit to 
the hadron-hadron data [0 yielded a universal pomeron trajectory ap{t) = ttjp(O) + a' t, with 
ajp{0) ~ 1.08 (the 'soft' pomeron). The dependence of the elastic hadron-hadron cross section 
on t is well described by an exponential distribution at small \t\, with a slope that increases 
with the cm. energy (shrinkage), leading to ~ 0.25 GeV~^. Regge theory can also be 
used to describe inclusive diffractive dissociation 0. These processes have been studied in 
hadron-hadron (see e.g. 0]-[0), photon-hadron ([T^-[T^) and DIS (0-0) interactions. 

For diffractive single dissociation in DIS, e^{k)p{P) e^{k')Xp{P'), the cross section can be 
expressed in terms of the diffractive structure function F.^^^^: 



l + {l-yy)Fr'{f3,Q\xp,t), (1) 



1 



where a is the electromagnetic couphng constant and the contributions of the longitudinal 
structure function and of Z° exchange have been neglected. Note that, with this definition, 
F^^^^ has dimensions of GeV~^. Integrating over t, we can define a structure function F^^^^: 



(fadiff 2Tra 
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The relevant kinematic variables are defined as follows. Deep inelastic scattering events are 
described by = —q^ = —{k — k'Y, the negative of the squared four-momentum transfer 
carried by the virtual photon; by the Bjorken variable x = -^^^^ and hy y = Q'^/xs, the 
fractional energy transferred to the proton in its rest frame, where y/s = 300 GeV is the ep 



cm. energy. The cm. energy of the virtual photon-proton (7*^) system is IV ~ yQ^{l/x — 1). 
Additional variables are required to describe diffractive scattering: 

t-{p-p'? X -i^^nn 3- - " (3) 

where x^ is the momentum fraction of the proton carried by the pomeron. In models where 



the pomeron has a partonic structure (see e.g. [0), /3 is the momentum fraction of the struck 
quark within the pomeron. 



Assuming factorization, as in the model of Ingelman and Schlein IT3, the structure function 
^2^^^^ is factorized into a pomeron flux, depending on x^ and t, and a pomeron structure 
function, which depends on (3 and Q^: 

F,''^'^ = Mx,,t)-Ff{l3,Q'). (4) 

In Regge theory the Xjp dependence of the flux can be expressed as (1/Xjp)^'^jp*^*^~^. 

In this paper, we present a measurement of the structure function F^^'^^ in the process ep — > eXp 
in the range 0.015 < P < 0.5, 4 x 10"'' < a;^, < 0.03 and averaged over between 5 and 20 
GeV^ and over t in the interval 0.073 < \t\ < 0.4 GeV^. We also present a measurement of the 
differential cross section as a function of t in a similar kinematic range. 



3 Experimental setup 
3.1 HERA 

In 1994 HERA operated with 153 colliding bunches of Ep = 820 GeV protons and E^ = 
27.5 GeV positrons. Additionally 15 unpaired positron and 17 unpaired proton bunches circu- 
lated and were used to determine beam related background. The integrated luminosity for the 
present study, which required the leading proton spectrometer to be in operating position (see 
below), is 900 ± 14 nb'^ 
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3.2 The ZEUS detector 



A detailed description of the ZEUS detector is given elsewhere [0]. The main components of 
the detector used in this analysis are the uranium-scintillator calorimeter, the tracking detectors 
and the leading proton spectrometer. 

The uranium-scintillator calorimeter (CAL) covers 99.7% of the solid angle and is divided 
into three parts, the forward^ (FCAL, covering the range 4.3 > rj > 1.1 in pseudorapidity) , 
the barrel (BCAL, 1.1 > t] > -0.75) and the rear (RCAL, -0.75 > t] > -3.8). Each part 
is longitudinally subdivided into electromagnetic (EMC) and hadronic (HAC) sections with 
typical cell sizes of 5 x 20 cm^ in the EMC (10 x 20 cm^ in the RCAL) and 20 x 20 cm^ in the 
HAC. The timing resolution is better than 1 ns for energy deposits greater than 4 GeV. The 



energy resolution was measured in test beams [19| to he a/E = 18%/vi?(GeV) for electrons 



and a/E = 35%/\'E{GeV) for hadrons. In order to minimize the effects of noise due to the 
uranium activity, the isolated EMC (HAC) cells with energy less than 100 (150) MeV were 
discarded from the analysis. 

The tracking system consists of two concentric cylindrical drift chambers, the vertex detector 
(VXD) [|20[ and the central tracking detector (CTD) [^, operating in a magnetic field of 
1.43 T. The CTD, which encloses the VXD, is a drift chamber consisting of 72 cylindrical 
layers, arranged in 9 superlayers. The measured resolution in transverse momentum for tracks 
with hits in all superlayers is <Jpj,/pT = 0.005pT 0.016 {pt in GeV). The interaction vertex is 
measured with a resolution of 0.4 cm in Z and 0.1 cm in the XY plane. 

The position of the scattered positron close to the rear beam pipe region is determined with 
precision by the small-angle rear tracking detector (SRTD), consisting of 2 planes of scintillator 
strips attached to the front face of the RCAL, covering an area of 68 x 68 cm^. A hole of 
20 X 20 cm^ at the centre of the RCAL and the SRTD accommodates the beam pipe. The 
SRTD signals resolve single minimum ionizing particles and provide a position resolution of 
0.3 cm. The time resolution is better than 2 ns for a minimum ionizing particle. The SRTD is 
also used as a presampler to correct the positron energy for losses in the inactive material in 
front of the rear calorimeter 



The proton remnant tagger (PRTl) [0 is a set of scintillator counters surrounding the beam 
pipe in the forward part of the ZEUS detector at Z = 5 m. The tagger consists of two layers of 
scintillating materials separated by a 1 mm thick lead absorber. Each layer is split vertically 
into two halves and each half is read out by a photomultiplier tube. The geometric acceptance 
of PRTl extends over the pseudorapidity region 4.2 < rj < 5.8. 

The luminosity is determined via the bremsstrahlung process e'^p —>■ e^p~f by measuring ener- 
getic photons in a lead-scintillator calorimeter (LUMI) placed aX Z = —107 m along the beam 
line [in. 



The LPS P3 measures protons scattered at very small angles with respect to the beam line and 



escaping the central apparatus through the forward beam hole. Such protons carry a substantial 



^The ZEUS coordinate system is right-handed with the Z-axis pointing in the proton beam direction, referred 
to as forward, and the X-axis horizontal pointing towards the centre of HERA. The pseudorapidity 7] is defined 
as — ln(tan |), where the polar angle 9 is measured with respect to the proton beam direction. 
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fraction Xl of the incoming proton momentum and have a small transverse momentum px 
(^1 GeV). The spectrometer consists of six detector stations located at distances of 23 to 
90 m along the proton beam line. In 1994 the three most forward stations S4, S5 and S6 were 
operational: each of these stations consists of an upper and a lower half, which partially overlap 
during data-taking. Each half is equipped with six rectangular parallel planes of silicon micro- 
strip detectors. Three different strip orientations (two vertical, two at +45° with respect to the 
vertical direction and two at —45°) are used, in order to remove reconstruction ambiguities. 
The dimensions of the planes vary from station to station and are approximately 4x6 cm^; the 
pitch is 115 /xm for the planes with vertical strips and 81 /xm for the other planes. The edges of 
the detectors closest to the beam have an elliptical contour which follows the lOcr profile of the 
beam, where a is the standard deviation of the beam spatial distribution in the transverse plane. 
During data-taking, the planes are inserted in the beam pipe by means of re-entrant Roman 
pots and are retracted during beam dump and filling operations of the HERA machine. The 
LPS coordinates are reconstructed with a precision of approximately 35 /im, which includes the 
intrinsic resolution of the coordinate reconstruction and the alignment precision of the detector 
stations. The track deflection induced by the magnets in the proton beam line is used for 
momentum analysis of the scattered proton. The xl resolution was determined using elastic 
p° photoproduction events to be better than 0.4% at 820 GeV; the intrinsic pt resolution is 
approximately 5 MeV at = 1 and is less than the spread introduced by the angular divergence 
of the proton beam (which is 40 MeV in the horizontal and 90 MeV in the vertical plane). 



4 Reconstruction of kinematic variables 

In order to reconstruct the kinematic variables x, and W (see section 2), the so-called 
"Double Angle" method was used, which derives the above quantities from the scattering 
angle of the positron and that of the struck (massless) quark. The latter is deduced from the 
momenta of all final state particles, except the scattered proton and positron. 

The final state particles in the reaction e^p — >■ e^Xp were reconstructed from the tracks and 
calorimeter energy deposits. The scattered positron identification algorithm was based on a 
neural network which uses information from the CAL. The momenta of the particles of 



the system X were reconstructed from calorimeter clusters and from tracks in the CTD. The 
calorimeter clusters were formed by grouping CAL cells into cones around the cell with a local 
energy maximum. The position of these objects was given by the sum of the positions of the 
single cells weighted by the logarithm of their energy. Reconstructed tracks were required to 
have transverse momenta of at least 100 MeV and were matched to clusters with a procedure 
based on the distance of closest approach. In case the cluster (track) was not matched to a 
track (cluster), the momentum of the particle was reconstructed only from the cluster (track). 
When the calorimeter cluster was matched to a track, the track momentum measurement was 
taken if the following two conditions were satisfied: the ratio between the energy of the cluster 
and the track momentum was less than 0.8; the track momentum resolution was better than the 
calorimeter energy resolution. If these two criteria were not satisfied, the particle momentum 
was reconstructed from the calorimeter cluster energy. An appropriate algorithm was developed 
for the cases in which more tracks pointed to a single cluster or more clusters were matched 
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to one track. The use of tracks improves the Mx resolution and reduces the sensitivity to the 
losses due to inactive material in front of the calorimeter. 



The quantity t was determined from the fractional momentum xl — p^^^ / Ep ~ p^^^/Ep of 
the scattered proton measured in the LPS and its transverse momentum p^^^ with respect to 
the proton beam direction: 

t=-<^. (5) 

XL 

The t resolution is approximately cr{\t\) = 140 MeVy^ {t in GeV^) and is dominated by the 
angular beam spread. 

The quantities and /3 were reconstructed using the values of and determined with 
the double angle method and the mass of the final hadronic system Mx- 

_ Ml + Q' 



The mass Mx was measured in two ways. In the first method, it was reconstructed as: 

{Mrn' = (E - {y.px? - i^P^r - {Y.pz)\ (7) 

h h h h 

where the sums run over all particles detected, except the scattered positron and proton, and 
P^ = (PxtPytPz) is momentum vector assigned to each particle of energy E^. From Monte 
Carlo studies the resolution on M^^"* was found to be approximately 40% for M^^"* < 3 GeV, 
18% in the range between 3 and 10 GeV and 14% for M^^°^ between 10 and 30 GeV. The 
mean reconstructed M^*'^"*' is 80% of its true value for all values of M™^"*: therefore an overall 
correction factor of 1.25 was applied. In the second method, Mx was reconstructed using the 
LPS and the double angle variables: 

{Mk^'^y = (1 - xl){W' + Q'- ml) -Q' + t, (8) 

where rUp is the proton mass. When xl approaches unity (low Mx), the resolution on M^^^ 
worsens. Monte Carlo studies showed a resolution of 17% for M^^^ > 10 GeV, of 35% for 
Mx^^ between 3 and 10 GeV and a substantially worse resolution for Mx^^ < 3 GeV. No 
additional correction factor was needed. 

In the analysis described here, the first reconstruction method was used to evaluate the mass 
Mx, calculated as 1.25 ■ M^'^"'^: the resulting resolutions on and f3 in the measured range 
are 30% and 25%, respectively. A mixed method with Mx^^ for high masses was used as a 
systematic check. 



5 Monte Carlo simulation and LPS acceptance 

The diffractive single dissociation process, e~^p e~^Xp, was modelled using two Monte Carlo 
(MC) generators which assume pomeron exchange, the RAPGAP program and a program based 
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on the Nikolaev and Zakharov model, as well as a MC generator for the exclusive reaction 



The RAPGAP p8| program is based on a factorizable model |]T^ in which the incoming proton 



emits a pomeron, whose constituents take part in the hard scattering. For the pomeron flux, 



the Streng parametrization |^ was used, which has an effective t dependence ranging from 
e~^'*' to e~^'*' in the range covered by the present measurement; for the pomeron structure 
function the form fitted in our previous analysis was taken, which contains a mixture of a 
"hard" (oc — (3)) and a "soft" (oc (1 — quark parton density, with no evolution. 
The program is interfaced to HERACLES |30| for the QED radiative corrections, the parton 



shower is simulated using the colour-dipole model as implemented in ARIADNE and the 
fragmentation is carried out with the Lund string model as in JETSET |32| . The region of low 
masses was tuned to the measured ratio of to p° resonance production |^. The events were 
generated with a lower xl cut of 0.9 and with |t| < 1 GeV^. 

The Nikolaev and Zakharov model pictures the diffractive dissociation of the photon as a 



fluctuation of the photon into a qq or qqg state |Q. The interaction with the proton proceeds 



via the exchange of a two-gluon colour-singlet state. The cross section can be approximated in 
terms of a two-component structure function of the pomeron, each with its own flux factor. The 
hard component, corresponding to the qq state, has a (3 dependence of the form /3(1 — /?) and 
an exponential t distribution with a slope of approximately 10 GeV~^. The soft component, 
which corresponds to the qqg state, has a /3 dependence of the form (1 — and a t slope of 
about 6 GeV~^. In the Monte Carlo implementation of this model [Q, hereafter referred to 
as NZ, the mass spectrum contains both components but the qqg states are fragmented (using 
JETSET) into hadrons as if they were a qq system with the same mass. The generation is 
limited to > 1.7 GeV. 

To improve the description of the low mass region, events in the region Mx < 1-7 GeV were 
generated separately using a MC for exclusive electroproduction, ep — > ep^p. In this MC 
the total cross section for the process 7*p —>■ p^p and the ratio of the longitudinal to the 
transverse cross section are derived from a parametrization of ZEUS and low energy data. 
The generator is interfaced to HERACLES which simulates the initial and final state radiation. 
The RAPGAP sample and the p° MC sample were mixed in such a way that the weighted sum 
reproduced the observed Mx distribution in the data: only a small contribution from the p° 
MC needed to be added and this mixture was used for the acceptance correction. In a similar 
way, the NZ and the p° MC were mixed and used as a systematic check. 

In order to estimate the background, the following processes were generated. The events 
where the proton dissociates into a state of mass Mjy {ep eXN) were generated with 
the EPSOFT |^ and PYTHIA programs, where the mass spectrum of the nucleon sys- 



tem and the ratio of double to single dissociation are generated according to measurements 
from proton-proton colliders f^. The spectrum is of the form 1/M^, with a ~ 2. To 
evaluate the background due to the one-pion exchange process, the RAPGAP MC implemen- 
tation of this exchange was used, in which the flux of pions is assumed to be fn{x-„,t) ~ 
(1 — 4m^)^/(l — 4t)^ ■ {xT,t/(t — m'lY) where m,r denotes the pion mass and x-„ is the 
fraction of the proton momentum carried by the pion. 
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All generated events were passed through the standard ZEUS detector simulation, based on the 
GEANT program , and through the trigger simulation and the event reconstruction package. 



For the scattered proton the simulation includes the geometry of the beam pipe apertures, the 
HERA magnets and their fields, and the response and noise of the LPS detectors. An accurate 
simulation of the interaction vertex position and of the effect of the proton beam tilt and 
emittance is also included. 

The distribution of the detected proton transverse momentum, p^^^, is limited by the ge- 
ometric acceptance of the LPS, which depends on the geometry of the beam pipe and the 
elliptical cutouts in the silicon detector planes. Figure 1 shows the LPS geometric accep- 
tance for tracks with xl = 1 in the px^^,pY^^ plane, where p^^ ^py^^ are the X and 
Y components of the proton momentum. Protons with GeV and ~ 1 are 

too close to the beam to be safely measured. In the region covered by the LPS planes 
(200 MeV< <400 MeV, |py^'^|;S600 MeV) the acceptance is large, as shown by the fines 

in tfie figure: however, integrating over the azimuthal angle (f) the resulting acceptance is ~ 6%. 
For Xl > 0.97 and 0.25<,P'^^'^<,0.65 GeV, the acceptance does not vary strongly with xl and 
changes from 15% at low p^^^ to 2% at high pff^^. The alignment uncertainty does not signif- 
icantly affect the geometric acceptance, and, in this analysis, cuts were applied to ensure that 
the events lie in a well understood region. 



6 Data selection 

Deep inelastic events in the ZEUS detector were first selected online by a three-level trigger 
system (details can be found in [M). At the first level DIS events were selected by requiring 
a minimum energy deposit in the electromagnetic section of the GAL. The threshold varied 
between 3.4 and 4.8 GeV depending on the position in the GAL. At the second level trigger, 
beam-gas background was reduced by using the measured times of energy deposits and the 
summed energies from the calorimeter. Events were accepted if S = J2i Ei{l—cos 9i) > 24 GeV— 
2E^, where Ei, 6i are the energies and polar angles of the calorimeter cells and E^ is the energy 
measured in the luminosity photon calorimeter, thereby accounting for the photon emitted in 
events with initial state radiation. For fully contained events, the quantity 5 is expected to be 
twice the positron beam energy, 5 ~ 55 GeV. At the third level, algorithms to reject beam-halo 
and cosmic muons, as well as a stricter 5 cut (5 > 25 GeV — 2i?^) , were applied, together with 
the requirement of a scattered positron candidate with energy greater than 4 GeV. 

Diffractive DIS candidate events were further selected offline in two steps: first the standard 
inclusive DIS selection was applied, then a high xl track in the LPS was required. 

Neutral current DIS events were selected as follows: 

• A scattered positron with E'^ > 10 GeV was required, where E'^ is the energy after 
presampler correction with the SRTD (when available). 

• The impact point of the scattered positron in the SRTD was required to be outside a 
square of 26 x 26 cm^ centred on the beam. This cut ensures full containment of the 
positron shower in RGAL. 
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40 < 5 < 65 GeV, to reduce the photoproduction background and the radiative correc- 
tions. 

The Z coordinate of the reconstructed vertex was required to be in the range —50 < Z < 
100 cm. 

T/e < 0.95, where y^, = 1 — E'^{1 — cos 9e)/2Ee is the value of y calculated from the positron 
variables {9e is the polar angle of the scattered positron). 

yjB > 0.03, where yjB = Y^h^hi^ — cos9h)/2Ee is the value of y calculated from the 
hadronic energy flow in this case the combination of tracks and clusters described 



in section 4 was used. 
The following cuts were applied to select diffractive events: 

• A track in the LPS was required, with 0.97 < xl < 1.01, where the lower limit was 
applied to reduce non-pomeron exchange contributions and to select a region of uniform 
acceptance, while the upper limit corresponds to a +2.5a distance from the xl = 1 peak, 
where a is the average LPS resolution at ~ 1. 

• The LPS track was extrapolated along the proton beam line. No track was accepted if, 
at any point, the minimum distance of approach to the beam pipe, Apip^, was less than 
500 fim. This cut reduces the sensitivity of the acceptance to the uncertainty in the 
position of the beam pipe apertures. 

• The total E + pz ^ {E + pz)ci^i + 2pY^ = Ei + cos 6*^) + 2pY^ of the event (the sum 
of the energy and the longitudinal component of the total momentum measured in the 
calorimeter and in the LPS) was required to be < 1655 GeV. For fully contained events 
this quantity should be equal to 2Ep = 1640 GeV; the cut, which takes into account the 
resolution on the measured value of Pz^^ , reduces the background due to overlay events 
(see following section). 

After this selection 553 events were left of which 376 were in the region between 5 and 
20 GeV^. 

Figures 2 to 4 show the distributions of the selected events and of the Monte Carlo model 
used in the analysis (weighted sum of the RAPGAP and the p° MC samples), as a function 
of Q'^,x,W,yjB (figure 2), of r]max, Mx,Xjp, P (figure 3) and of XL,\t\,Px^^ ,Py^^ (figure 4). 
The distributions are all uncorrected and the MC histograms are normalized to the number of 
events in the data; a cut of > 5 GeV^ was apphed to both data and Monte Carlo samples. 
In fig. 3a, the variable rj^ax is the maximum pseudorapidity of all calorimeter clusters with an 
energy of at least 400 MeV or tracks with momentum of at least 400 MeV. In our previous 
analysis, events with a large rapidity gap were defined by rimax^^-5-2.5 [Q: however the events 
with a tagged proton with xl > 0.97 can have larger values of rjmax- This allows us to study 
events with Mx^20 GeV, where the hadronic system extends close to the proton beam pipe, 
and no gap may be observed in the detector. In figs. 3c and 3d, the In distribution is shown 
in two W bins {50 < W < 120 GeV, 120 < W < 270 GeV). The data and the diffractive 
MC model are in good agreement. Compared to our previous analyses, the data extend the 
explored kinematic region to higher values of Mx, and therefore to higher Xp and lower (3. 
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7 Background 



After the selection described in the previous section, the sample still contains some background, 
mainly due to beam-halo and DIS processes which are not single diffractive dissociation. 

Proton beam-halo events originate from interactions of beam protons with the residual gas in 
the pipe or with the beam collimators. These events have a scattered proton of energy close to 
that of the beam: when they accidentally overlap with a genuine DIS event, they may give a 
false diffractive signal. Most of these events, however, appear to violate energy and momentum 
conservation. Figure 5 shows a scatter plot of xl and the {E + pz)cai (measured with the 
calorimeter) for the DIS events selected in the LPS (excluding the xl and the E + pz cuts): 
a clear band at xl ~ 1, uncorrelated with the energy measured in the calorimeter, can be 
ascribed to beam-halo events. This type of background can be rejected by requiring that the 
total E + Pz of the event be conserved, as for beam-halo it can exceed the kinematic limit of 
1640 GeV (see line in fig. 5). To evaluate the residual background after the cuts mentioned 
in the previous section, the 2p^^^ distribution of events with unphysical tracks in the LPS 
{{E + pz)ca.\ + 2p^^^ > 1655 GeV) was randomly mixed with the {E + pz)cai distribution 
for DIS events, to create a. E + pz distribution for beam-halo events. The obtained E + pz 
distribution was normalized to the observed data distribution for E + pz > 1655 GeV. The 
remaining background (below the E + pz = 1655 GeV cut) was estimated to be less than 6%. 

In order to evaluate the background due to proton dissociation and pion exchange, the cut on 
Xl was removed. Figure 6 shows the observed uncorrected Xl spectrum for the data (dots): 
a narrow diffractive peak is seen at xl — 1, together with a distribution at lower xl due to 
the background processes mentioned above. Note that the acceptance falls by almost an order 
of magnitude between xl — 0.8 and xl — 1- These processes were modelled using the Monte 
Carlo programs described in section 5, EPSOFT for the proton dissociation, RAPGAP for 
the pion exchange, while for the diffractive signal the combined RAPGAP and MCs were 
used. The weighted sum of the three components was fitted to the observed xl spectrum in 
the data in two steps. The weight factor for the proton dissociation events was first determined 
by normalizing the EPSOFT MC sample to the data in the region xl < 0.95 and rjmax < 1-5, 
where double dissociation dominates. The weight factors for the single diffraction and the pion 
exchange processes were determined by fitting the weighted sum of the three MC samples to the 
observed xl distribution. The resulting sum is shown as the solid line in fig. 6. After applying 
the Xl > 0.97 cut, the background due to proton dissociation is less than 3%, to be compared 
with the 10-15% contamination estimated to be present in our previous analyses The 
fit was repeated using different MC models (PYTHIA for proton dissociation, NZ plus p° for 
the signal) and consistent results were found. This estimate of the background due to proton 
dissociation was checked by looking at events tagged both in the proton remnant tagger PRTl 
and in the LPS. According to simulations based on the EPSOFT Monte Carlo, the PRTl has a 
tagging efficiency of approximately 50% for proton dissociation events; this leads, after the xl 
cut, to an estimated contamination below 3%, consistent with the evaluation described above. 

The background due to the pion exchange was found to be less than 1% after the xl > 0.97 cut. 
Note that the fit to the xl spectrum is used to give an estimate of the upper limit on the pion 
exchange background at high xl and it is not meant to be a complete study of the distribution. 
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The background due to other reggeon exchanges was not included in this evaluation: at high 
Xjp values (x^p >;0.01), a contribution from other reggeon exchanges is likely however the 
predictions vary significantly in different models E^, The possible contribution of these 



additional trajectories at high is discussed in section 9. Background due to non-diffractive 
DIS processes (not shown in the figure), in which one of the proton fragments is observed in 
the LPS, gives a contribution at small values of xl- 

The background due to non-ep interactions (excluding beam-halo events) was evaluated to be 
negligible from the data taken with the unpaired proton bunches. The background due to 
photoproduction events was found to be negligible from Monte Carlo studies. 

As the backgrounds due to beam-halo, non-ep interactions, proton dissociation and pion ex- 
change were found to be small compared with the statistical precision of the data, they were 
not subtracted in the results shown in the following. 



8 Measurement of the t distribution 

The measurement of the t dependence in diffractive DIS was limited to the kinematic range 
0.073 < \t\ < 0.4 GeV2, 5 < < 20 GeV^, 50 < W < 270 GeV, 0.015 < /3 < 0.5 and 
Xl > 0.97. The range in \t\ was limited to values greater than 0.073 GeV^, since for lower 
values of \t\ the acceptance varies rapidly; the upper limit of 0.4 GeV^ restricts the data to a 
region where the LPS acceptance exceeds 2%. The bin widths in t were chosen to be larger than 
the resolution, resulting in four bins. The acceptance and the detector effects were unfolded 
using a bin-by-bin correction determined with the RAPGAP plus p° MC. 

Figure 7 shows the measured differential cross section, da/dt. The distribution was fitted with 
a single exponential form, shown as the solid line: 

da/dt = Ae^'; (9) 

the value of the fitted parameter b is: 

b = 7.2± l.l(stat.)+[;;9(syst.) GeV'^, (10) 

where the first uncertainty is statistical and the second is systematic. 

The systematic uncertainties can be subdivided into three groups, those due to the DIS se- 
lection, those due to the LPS acceptance and background and, finally, those related to the 
unfolding. The systematic error due to the DIS acceptance was evaluated by changing the fol- 
lowing cuts: the positron energy cut was moved to 8 GeV and 12 GeV, the cut on the impact 
position of the positron using the SRTD was moved to 24 x 24 cm^ and to 28 x 28 cm^; the yjB 
cut was changed to 0.02 and 0.04. The effects on the slope b varied between +5.5% and —3%. 
To estimate the systematic contributions due to the LPS acceptance, the following checks were 
performed: the Apjpe cut was increased to 0.1 cm causing a negligible variation to the slope. 
The values of pjf^ and py^^ were shifted by ib3 MeV and ±6 MeV, respectively, as a test of 
the influence of the alignment procedure (see ||2^), yielding variations on b between —6% and 
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+3%. The mass Mx was reconstructed using the calorimeter for low masses and the LPS for 
high masses, giving changes around ±5% on the slope. The NZ plus p° Monte Carlo was used 
for the acceptance calculation, the t slope in RAPGAP was changed by ±1 GeV~^, the a'^ 
value in the MC was changed from 0.25 GeV~^ to 0. These changes lead to variations of at 
most 4% for the slope h. The total systematic error was obtained by adding in quadrature the 
positive and negative deviations separately. 

In a Regge-type approach P], the slope of the exponential t distribution in single diffractive 
interactions is predicted to be 6 ^ 6o + 2a;^ ln(l/x^). Assuming the value of 6o = 4.6 GeV~^ in- 
ferred by Goulianos from elastic pp scattering at y/s = 1800 GeV |^ and a'^ = 0.25 GeV~^ 
10, in the present kinematic range (< >^ 8 x 10"'^) the measured b value is compatible 
with the predicted value of ~ 7 GeV~^ for a soft pomeron exchange; it is also consistent with 



the values predicted by some perturbative QCD models |47 . 



9 Measurement of F^^^^ and F^^^^ 

The measurement of the proton in the LPS permits the determination of the cross section for 
single diffractive dissociation in DIS as a function of the four kinematic variables j3, Q^, and t. 
Given the small statistics collected in the 1994 running, the measurement was performed in a 
single t bin, 0.073 < \t\ < 0.4 GeV^ and a single Q"^ bin, 5 < < 20 GeV^, with average 
values of < \t\ >= 0.17 GeV^ and < >= 8 GeV^. The sizes of the (3 and x^ bins were 
chosen to be larger than the resolutions (see section 4); a minimum number of 8 events in 
each bin was also required. The chosen bins are shown in figure 8: they are in the range 
4 X 10~^ < Xp < 0.03, 0.015 < P < 0.5. The acceptance, which includes the efficiency of 
the DIS selection and the geometric acceptance of the LPS, varies between 3% and 11% in 
these bins. The purity, defined as the ratio of the number of MC events generated in a bin 
and reconstructed in the same bin over the number of events reconstructed in that bin, was 
required to be greater than 25% in the bins used for the measurement and is typically more 
than 40%. The effect of the longitudinal structure function Fl is assumed to be smaller than 
the statistical errors in the kinematic range considered ||^ and is neglected. 

In order to extract the structure function F^^^\ the weighted sum of the RAPGAP and the p° 
MC, which gives a good description of the data, was used to obtain the correction for acceptance 
with a bin-by-bin unfolding method. In addition a bin centring correction was applied, using 
the F^^"^^ parametrization of the RAPGAP MC. 

The following systematic checks were performed for the measurement of F^^^^: 

• To check the photoproduction background and the QED radiative corrections, the cut on 
6 was changed to 37 GeV and 42 GeV, and the cut on the scattered positron energy was 
varied to 8 GeV and 12 GeV. The effect was typically up to 3%. 

• To check the acceptance at low which is determined by the positron position, the cut 
on the impact point of the positron was changed to 24 x 24 cm^ and to 28 x 28 cm^; both 
induced variations on F^ smaller than 15%. 
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An alternative method to reconstruct the kinematics, the S method [|8|, was used, giving 
typical variations of about 20%. 

The cut on yjs was moved by ± 0.01, yielding changes of typically around 10%. 

The cut on the minimum energy deposit of the EMC (HAC) cells was moved to 140 
(160) MeV and to 80 (120) MeV, to check the effect of the calorimeter noise on the mass 
reconstruction. The effect was between 10% and 25% for the first cut, and up to 15% for 
the second check. 

The LPS track selection was modified by raising the Apjpe cut to 0.1 cm, producing 
variations of F^^^^ of less than 12%. 



To check the background estimation due to one-pion exchange and proton dissociation, 
MC predictions were statistically subtracted from the bins, producing changes up to 5% 
in a few bins. The beam-halo was also subtracted with negligible effects on F2 ■ The 
xl cut was moved to 0.96, causing changes on F2^'^'' of less than 5%, except in one bin 
at high Xjp where the change was 15%. 

An alternative method to reconstruct the mass Mx was employed, which uses 1.25 -M^p*^"* 
at low masses and Mx^^ at high masses, giving an effect of less than 2% at high /3 and 
up to 25% at low [3. 



An unfolding method based on Bayes' theorem [|9| was applied for the acceptance cor- 
rection causing variations typically less than 3%. The NZ plus p° Monte Carlo was used 
to unfold, giving changes typically less than 20%, except at low j3 where the changes 
were up to 40%. The Xjp and t dependence in RAPGAP were reweighted to a fixed t 
dependence of the form (l/a;jp)e~^l*l or to a fixed x^ dependence of the form (l/x^)^'^, 
yielding changes of less than 5%. 



Most of the systematic checks yielded results which agree with the central value within the 
statistical errors. The negative and positive deviations in each bin were separately combined 
in quadrature. 

The results for F2^^'' are given in table 1. The systematic errors do not include a 5.5% overall 
normalization uncertainty due to the luminosity determination and trigger efficiency (2%) and 
due to the LPS acceptance varying because, due to the proton beam conditions, the LPS stations 
had to be positioned differently in different runs (5%). The values of ^2^*^^'' are plotted in fig. 9 
as Xjp ■F2^^^ as a function of x^, in four (3 intervals, with central values oi j3 = 0.028, 0.07, 0.175 
and 0.375, respectively. The ^2^^^"* data are observed to fall rapidly with increasing Xp. To 
investigate whether the x^ dependence is the same in all j3 intervals, we performed fits of the 
form: 

F2^(') = A, . il/x^Y , (11) 

where the normalization constants Ai were allowed to vary in each bin, while the exponent a 
was constrained to be the same in all four (3 bins. The result of the fit is 

a(0.073 < \t\ < 0.4 GeV^) = 1.00 ± 0.09 (stat.)+J;JJ(syst.), (12) 
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with a value of 10 for 8 degrees of freedom, showing that the result is consistent with a single 
Xjp dependence in all (3 bins. The systematic error was obtained by re-fitting the -Fj^'''^'* values 
obtained for each of the systematic checks mentioned above. The main contributions arose 
from changing the unfolding procedure (4%), using the S method (5%) and applying different 
noise suppression cuts (7%). 

In order to compare with previous measurements of the diffractive structure function, we ex- 
tracted F^^^\j3 , Q"^ , X jp) by integrating F^^^^ over t. The extrapolation to the whole t range 
(0 < \t\ < oo) was performed using the F^^^^ parametrization in the MC. The result is shown 
in table 1 and is plotted as ■ F^^^^ in figure 10. Fitting our data in each (3 interval with the 
form {l/XpY yielded a value a = 1.01 ± 0.10 (stat.)^g Qg(syst.). 

This value of a is much lower than the value obtained with the ZEUS 1993 data based 
on a large rapidity gap analysis^. The ZEUS 1993 F.^^^^ result corresponds to 0.1 < /5 < 0.8 
and 6.3 x 10""^ < x^, < 10~^, for 8 < Q"^ < 100 GeV^, and it is compatible with a single x^ 
dependence in all (3 bins, with slope a = 1.30 ± 0.08(stat.)l'iQ°4(syst.). As already mentioned, 
the present analysis covers a different kinematic range, extending to lower (3 and higher x^ at 
< >= 8 GeV^. The lower value of a compared to our previous result may be ascribed to 
the presence of additional subleading trajectories contributing in the x^ range covered by this 
analysis. This contribution has been already observed in the analysis of the HI 1994 data 
which are shown for comparison in fig. 10. Only the points at the and f3 values closest to the 
ones presented here are shown. Assuming pomeron exchange to be the dominant contribution, 
the values of a(0.073 < \t\ < 0.4 GeV^) and a can be related, respectively, to the value in 
the given t range and U^, the pomeron trajectory averaged over t (see section 2). However, 
as already mentioned in section 7, at high x^p a contribution from reggeon exchange cannot be 
excluded, in addition to pomeron exchange: the statistical precision of our data does not allow 
us to identify this component. As the reggeon trajectory has an intercept — 0.5, leading to 
an Xjp slope a ~ 0, its presence would lower the expected x^ slope in the x^ region where the 
reggeon exchange contributes. The ge covered by this analysis could then explain the 

different result on the slope d compared to our previous result. Moreover, models where the 
pomeron component is not factorizable can lead to a different x^ slope parameter, depending 
on the kinematic range. 

We have also examined the (3 dependence of the diffractive structure function. In order to do 
this, the F^^^^ values within each (3 bin are plotted in figure 11, extrapolated to a common 
value of Xjp = 0.01, assuming a universal, fixed slope a = 1.01. The solid line represents a fit of 
the form b ■ {1/Xjp)^[(3{l — (3) + c/2(l — (3^] to the measured F^^^\ The values obtained in the 
fit are: b = 0.087 ± 0.015 (stat.)+J;j2j(syst.) and c = 0.34 ± 0.11 (stat.)+|5;?o(syst.). This fit to 
the present data indicates that both a hard (oc l3{l — (3)) and a soft (oc (1-/9)^) component are 
needed. The parametrization obtained in our previous analysis 0|, extrapolated to Xjp = 0.01 
and scaled down by 15% to take into account the proton dissociation background contribution, 
is also shown in the figure as the dashed curve. 

^ We do not compare here to the result in ref . [|j , as this result was affected by a technical error in the Monte 
Carlo generation used for the unfolding, which led to a higher a value by about one unit of the quoted error. 
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10 Conclusions 



Diffractive DIS events have been studied at HERA using the ZEUS leading proton spectrometer. 
A clean sample of events was selected by requiring a scattered proton with xl > 0.97. The 
background due to proton dissociation was estimated to be approximately 3%, substantially 
lower than in our previous analyses. The use of the LPS has also allowed an extension of the 
ZEUS measurements to values of the final state hadronic mass Mx as high as 35 GeV. 

The t dependence was measured for the first time in this process in the range 0.073 < \t\ < 
0.4 GeV^, 5 < < 20 GeV^, 0.015 < /3 < 0.5 and 50 < < 270 GeV. The resulting 
distribution is described by the function e***, with b — 7.2± l.l(stat.)^Q g(syst.) GeV"^. 

The diffractive structure function F^^'^\p, Q^, Xjp,t) was measured in the interval 0.015 < (5 < 
0.5, 4 X 10~^ < < 3 X 10~^ and averaged over the range 0.073 < \t\ < 0.4 GeV^ and 
5 < < 20 GeV^. Because of the hmited statistical precision of the data, it is not possible 

to determine whether a different dependence is needed in different x^ and f3 ranges. The 
Xjp dependence is consistent, in all f3 intervals, with the form (l/x^)", with a(0.073 < \t\ < 
0.4 GeV^) = 1.00 ± 0.09 (stat.)^|]'g5(syst.). Integrating over t, the structure function F^^^^ was 
determined. A fit of the form (1/x^)" to ^2^^^^ yielded a = 1.01 ± 0.10 (stat.)+U;Jg(syst.). 

The result for the effective .x^ slope is lower than that obtained in our previous measurement. 
This analysis, however, extends the ^2^*'^'' measurement to values of Xp up to 0.03, where a 
significant component of reggeon exchange could contribute to lowering the effective slope 
parameter a. 



Acknowledgements 

We thank the DESY Directorate for their strong support and encouragement. We acknowledge 
the assistance of the DESY computing and networking staff. We are very grateful to the 
HERA machine group: collaboration with them was crucial to the successful installation and 
operation of the LPS. We would like to thank B. Hubbard for his invaluable contributions to 
the experiment, and the LPS in particular. It is also a pleasure to thank N.N. Nikolaev and 
M.G. Ryskin for useful discussions. 

References 

[1] ZEUS Collaboration, M. Derrick et al., Phys. Lett. B315 (1993) 481. 
[2] HI Collaboration, T. Ahmed et al., Nucl. Phys. B429 (1994) 477. 
[3] HI Collaboration, T. Ahmed et al., Phys. Lett. B348 (1995) 681. 
[4] ZEUS Collaboration, M. Derrick et al., Z. Phys. C68 (1995) 569. 

14 



[5] HI Collaboration, C. Adloff et al., DESY report DESY 97-158. 

[6] ZEUS Collaboration, M. Derrick et al., Z. Pliys. C70 (1996) 391. 

[7] A. Donnachie and P.V. Landshoff, Nucl. Phys. B244 (1984) 322; 
Phys. Lett. B296 (1992) 227. 

[8] A.H. Mueller, Phys. Rev. D2 (1970) 2963; ibid. D4 (1971) 150; 
D. Field and G. Fox, Nucl. Phys. B80 (1974) 367. 

[9] For a review, see e.g. K. Gouhanos, Phys. Rep. 101 (1983) 169. 

[10] UA4 Collaboration, D. Bernard et al., Phys. Lett. B186 (1987) 227. 

[11] UA8 Collaboration, R. Bonino et al., Phys. Lett. B211 (1988) 239; 
UA8 Collaboration, A. Brandt et al., Phys. Lett. B297 (1992) 417. 

[12] E710 Collaboration, N.A. Amos et al, Phys. Lett. B301 (1993) 313. 

[13] CDF Collaboration, F. Abe et al., Phys. Rev. D50 (1994) 5535. 

[14] E612 Collaboration, T.J. Chapin et al., Phys. Rev. D31 (1985) 17. 

[15] HI Collaboration, C. Adloff et al., Z. Phys. C74 (1997) 221. 

[16] ZEUS Collaboration, J. Breitweg et al., Z. Phys. C75 (1997) 421. 

[17] G. Ingelman and P. Schlein, Phys. Lett. B152 (1985) 256. 

[18] ZEUS Collaboration, M. Derrick et al, Phys. Lett. B293 (1992) 465; 
The ZEUS Detector, Status Report 1993, DESY 1993. 

[19] A. Andresen et al, Nucl. Inst. Meth. A309 (1991) 101; 
A. Caldwell et al., Nucl. Inst. Meth. A321 (1993) 356; 

A. Bernstein et al., Nucl. Inst. Meth. A336 (1993) 23. 

[20] C. Alvisi et al., Nucl. Inst. Meth. A305 (1991) 30. 

[21] N. Harnew et al., Nucl. Inst. Meth. A279 (1989) 290; 

B. Foster et al., Nucl. Phys. B (Proc. Suppl.) 32 (1993) 181; 
B. Foster et al., Nucl. Inst. Meth. A338 (1994) 254. 

[22] ZEUS Collaboration, M. Derrick et al., Z. Phys. C69 (1996) 607. 

[23] J. Andruszkow et al., DESY Report DESY 92-066; 

ZEUS Collaboration, M. Derrick et al., Z. Phys. C63 (1994) 391. 

[24] ZEUS Collaboration, M. Derrick et al, Z. Phys. C73 (1997) 253. 

[25] U. Amaldi et al., Phys. Lett. B43 (1973) 231. 



15 



[26] S. Bentvelsen, J. Engelen and P. Kooijman, Proceedings of the Workshop "Physics at 
HERA", vol. 1, DESY (1992) 23. 



[27; 
[28 
[29 

[30 

[31 

[32 
[33 
[34 

[35 
[36 
[37; 
[38 

[39 
[40 
[41 

[42' 

[43 



[44; 



[45 
[46 
[47 



A. Abramowicz, A. Caldwell and R. Sinkus, Nucl. Inst. Meth. A365 (1995) 508. 
H. Jung, DESY Report DESY 93-182. 

K.H. Streng, in Proc. of the Workshop "Physics at HERA", p. 365, ed. R.D. Peccei (Ham- 
burg 1987). 

K. Kwiatkowski, H. Spiesberger and H.-J. Mohring, Proceedings of the Workshop "Physics 
at HERA" vol. 3, DESY (1992) 1294. 

L. Lonnblad, Comp. Phys. Comm. 71 (1992) 15; 
L. Lonnblad, Z. Phys. C65 (1995) 285. 

T. Sjostrand and M. Bergtsson, Comp. Phys. Comm. 43 (1987) 367. 

ZEUS Collaboration, M. Derrick et al, Phys. Lett. B380 (1996) 220. 

N.N. Nikolaev and B.C. Zakharov, Z. Phys. C53 (1992) 331; 

M. Genovese, N.N. Nikolaev and B.C. Zakharov, JETP 81 (1995) 625. 

A. Solano, Ph.D. Thesis, University of Torino 1993 (unpublished). 

ZEUS Collaboration, M. Derrick et al, Phys. Lett. B356 (1995) 601. 

M. Kasprzak, PhD thesis, Warsaw University (1996), DESY Internal Report F35D-96-16. 

H.-U. Bengtsson and T. Sjostrand, Comp. Phys. Comm. 46 (1987) 43; 
T. Sjostrand, CERN TH-7112-93, (1994). 

L.L. Frankfurt, L. Mankiewics and M.I. Strikman, Z. Phys. A334 (1989) 343. 
GEANT 3.13, R. Brun et al., CERN DD/EE/84-1 (1987). 
ZEUS Collaboration, M. Derrick et al, Z. Phys. C72 (1996) 399. 

F. Jacquet and A. Blondel, Proceedings of "The study of an ep facility for Europe" , DESY 
79/48 (1979) 391. 



K. Golec-Biernat, J. Kwiecihski and A. Szczurek, INP Cracow preprint 1754/PH, |hep- 



ph/9701254 . 



N.N. Nikolaev, W. Schafer and B.C. Zakharov, preprint KFA-IKP(Th)-1996-06, [Eep^ 
ph/9608338| 



K. Goulianos, Phys. Lett. B358 (1995) 379 (and Erratum in Phys. Lett. B363 (1995) 268). 
CDF Collaboration, F. Abe et al., Phys. Rev. D50 (1994) 5518. 
N.N. Nikolaev, B.C. Zakharov and V.R. Zoller, Phys. Lett. B366 (1996) 337. 

16 



[48] U. Bassler and G. Bernardi, Nucl. Inst. Meth. A361 (1995) 197. 
[49] G. D'Agostini, Nucl. Inst. Meth. A362 (1995) 487. 



17 





Xjp 


Nobs 


F2^ ' ± Stat. ± sys. GeV^^ 


r^D(3) , , , , 

F^^ ' ± Stat. ± sys. 


0.028 


0.011 


14 


2.8 ± 0.8^^:^ 


1.3 ± 0.4+^:;^ 


0.028 


0.024 


17 


2.2 ± 0.6tii 


1.0 ± Q.stli 


0.07 


0.0044 


13 


7.8 ± 2.3^1? 


3.6± l.ltl;|] 


0.07 


0.011 


22 


3.8 ± 0.9l{i 


1.8 ± 0.4;[|;^ 


0.07 


0.024 


13 


1.8 ± 0.5l°i 


0.8±0.2^[j:f 


0.175 


0.0018 


15 


28.4 ± 7.71^8^ 


13.4±3.6t^i 


0.175 


0.0044 


21 


10.8 ±2.5^?:^ 


b.O±1.2t'oi 


0.175 


0.011 


22 


5.3 ± 1.211:1 


2.5 ± 0.6t|j| 


0.175 


0.024 


13 


3.6± i.itl:? 


1.7 ± 0.5lg:? 


0.375 


0.0007 


19 


89.1 ± 21.7t\'li 


42.7± 10.4l|^ 


0.375 
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53.6 ± 10.4l?i62 


25.3 ±4.91^:^ 


0.375 


0.0044 


16 


10.4 ± 2.7^^ 


4.9 ± i.3t;i 


0.375 


0.011 


8 


3.9 ± lAtli 


1.8±0.7tj;° 



Table 1: Results on Fp^> and F2 ^'^^ using the ZEUS LPS 1994 data at = g GeV^. The 
Fo values are in the t interval 0.073 < \t\ < 0.4 GeV^, the ^2^^^^ values are in the t range 
< |i| < 00. The overall normalization uncertainty of 5.5% is not included and no background 
is subtracted from the data. 
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Figure 1: Geometric LPS acceptance for tracks with xl 
indicate the regions where the acceptance is greater than 5%, 50% and 95%. 



1 in the Px^'^ ,Py^^ plane. The hnes 
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ZEUS 1994 




Figure 2: Observed event distributions as a function of (a) Q^, (b) x, (c) W, and (d) yjs of the 
reconstructed data (dots) compared to the Monte Carlo model (solid line). The Monte Carlo 
is the weighted sum of the RAPGAP and p° samples. The errors, shown as vertical bars, are 
statistical only. 
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ZEUS 1994 




Figure 3: Observed event distributions as a function of (a) rjmax , (b) Mx, (c) lnM|^ for 
50<W < 120 GeV, (d) InM^ for 120 < 1^ < 270 GeV, (e) and (f) P of the reconstructed 
data (dots) compared to the Monte Carlo (sohd hne). The Monte Carlo is the weighted sum 
of the RAPGAP and p° samples. The errors are statistical only. 
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Figure 4: Observed event distributions as a function of (a) Xi^ , (b) \t\, (c) pjf^ and (d) py^^ 
of the reconstructed data (dots) compared to the Monte Carlo (sohd hue). The Monte Carlo 
is the weighted sum of the RAPGAP and p° samples. The errors are statistical only. 
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Figure 5: A scatter plot of a;^ and the {E + pz)cai as measured in the calorimeter in the 
data (before the xl and the E + Pz cut). The region below the solid hne indicates the allowed 
kinematic region. 
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ZEUS 1994 
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Figure 6: The observed xl spectrum in the data (dots) where the xl cut at 0.97 has been 
removed. Overlaid is the resuh of fitting the distribution with a sum (full line) of the contri- 
bution due to proton dissociation (EPSOFT MC - shaded area), of the maximum contribution 
due to pion exchange (dotted line) and of the single photon dissociation signal (RAPGAP plus 
p° MC - dashed line). 
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ZEUS 1994 




Figure 7: The differential cross section da/dt for diffractive DIS events with a leading proton 
of XL > 0.97, in the range 5 < < 20 GeV^, 50 < < 270 GeV and 0.015 < p < 0.5. 
The inner error bars indicate the statistical errors, the outer error bars show the statistical and 
systematic errors added in quadrature. The line is the result of the fit described in the text. 
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Figure 8: The distribution of the events used in the analysis in the /3, plane. The bins in 
(3, Xjp used for the extraction of ^2^''''^ and F^f are shown. 
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Figure 9: The structure function ■ F2 {(3,Q'^,Xjp,t), plotted as a function of Xjp in four 
f3 bins, at < >— 8 GeV^ and < \t\ >— 0.17 GeV^. The inner error bars indicate the 
statistical errors; the outer error bars are the sum of the statistical and systematic errors added 
in quadrature. The 5.5% normalization uncertainty is not included. The solid line corresponds 
to the fit described in the text. 
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Figure 10: The function x^^ ■ F2 {(3,Q'^,Xjp) for this analysis at < >= 8 GeV^ and 
< \t\ < 00, compared to the results at < >= 7.5 GeV^ and |t| < 1 GeV^ from ref. [5]. 
The statistical and systematic errors are added in quadrature. The 5.5% and 6% normalization 
uncertainties of the ZEUS 1994 LPS data and of the HI 1994 data, respectively, are not included. 
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ZEUS 1994 




Figure 11: The value of Fp^\/3, < >^ 8 GeV^, = 0.01), obtained from the fits to the 
individual bins in j3. The solid line represents the fit to the data as described in the text. The 
dashed line indicates the parameterization of ref. [4] scaled down to remove the estimated 15% 
double dissociative contributions. 
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